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Preface 


PEEK  will  fly  in  an  F-18  Hornet  this  year.  Is  PEEK 
really  that  familiar  now?  The  answer  is  yes  and  no.  This 
research  bridged  the  gap  between  the  more  commonly  under¬ 
stood  mechanical  properties  of  the  PEEK  matrix  and  graphite 
fiber  composite  called  APC-2  and  the  less  well  understood 
viscoelastic  characteristics  that  it  exhibits. 

The  flexure  model  developed  here  is  a  simple  one,  with 
viscoelastic  parameters  substituted  for  elastic  material 
constants  under  assumed  quasi-static  conditions.  Although 
the  results  are  encouraging,  my  test  technique  can  be  im¬ 
proved,  and  these  recommendations  are  respectfully  passed  on 
to  the  Materials  Laboratory. 

Without  the  sponsorship  of  the  Non-Metallics  Division, 
Structural  Materials  Branch,  this  study  wouldn't  have  been 
possible.  I  am  indebted  to  my  sponsor.  Dr.  Charles  Lee,  and 
to  his  colleagues.  I  would  like  to  single  out  the  support 
of  Mr.  Ken  Lindsay,  Mr.  Bill  Ragland,  and  Ms.  Margie  Hudock . 

Without  the  help  of  the  faculty  at  APIT,  I  would  have 

misunderstood  the  importance  of  careful  analysis.  Special 

thanks  to  my  advisor.  Major  Paul  Copp,  for  the  concern  that 

he  showed  throughout  this  year.  _ 

?sslon  For 

Finally,  this  research  would  have  been  joyless  had  it  gram  — 

’  TAB 

not  been  for  Paige,  my  wife.  You  kept  the  big  picture  near 

1  an. 

me  at  all  times. 

Richard  P.  Lockwood 
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Abstract 


'The  thermoplastic  polyetheretherketone  (PEEK)  is  being 
considered  for  structural  applications  in  the  Advanced  Tac¬ 
tical  Fighter.  Previous  efforts  characterizing  PEEK  mechan¬ 
ical  properties  are  reviewed.  A  study  developing  the  visco¬ 
elastic  f requency-temperature  response  to  three  point 
flexure  of  neat  resin  and  PEEK/graphite  fiber  composites  is 
detailed.  The  behavior  of  unidirecti onal  laminates  is  com¬ 
pared  to  predictions  from  elastic  plate  theory,  which  also 
combines  micro-  and  macromechanical  approaches  to  predict 
transverse  and  longitudinal  complex  moduli.  The  experimen¬ 
tal  frequency-temperature  shift  factors,  obtained  under  low 
strain  in  a  Dynastat  dynamic  viscoelastic  analyzer,  are 
compared  to  shift  factors  predicted  from  the  Williams- 
Landel-Ferry  equation.  Determining  composite  behavior  from 
knowledge  of  matrix  and  fiber  behavior  is  still  applicable. 
The  capability  of  the  Dynastat  to  perform  viscoelastic 
characterizations  of  stiff,  unidirectional  composites,  based 
on  comparisons  of  the  various  shift  factors,  is  reduced  due 
to  the  presence  of  the  reinforcing  fibers.  Recommendations 


for  future  test  and  evaluations,  including  increasing  the 
dynamic  s  ignal-to-noise  ratio,  are  descr  ibed .  v  t  ./ 


VISCOELASTIC  BEHAVIOR  OF  A  POLYETHERETHERKETONE 

(PEEK)  COMPOSITE 

I .  Introduction 

Advances  in  aerospace  vehicle  performance  often  result 
from  the  improvements  in  the  mechanical  characteristics  of 
the  new  materials  in  the  vehicle.  The  design  process  re¬ 
quires  material  parameters  in  the  stress  analysis  and  in  the 
study,  of  failure  mechanisms.  This  research  characterizes 
the  viscoelastic  behavior  of  a  thermoplastic  and  graphite 
fiber  composite.  The  applicability  of  the  three-point 
flexure  to  this  characterization  is  considered,  along  with 
the  effects  of  viscoelasticity  on  assumptions  made  in  the 
development  of  an  analysis  of  composite  flexure. 

The  recent  development  of  polyetheretherketone  (PEEK), 
and  its  consideration  for  aerospace  structural  applications, 
results  from  the  improved  strength,  toughness,  stiffness, 
ease  of  producibil i ty,  resistance  to  solvents,  and,  of 
course,  the  light  weight  exhibited  by  graphite  reinforced 
composites  of  the  polymer.  Some  of  the  mechanical  proper¬ 
ties  already  known  about  this  material  are  described  below. 
Because  PEEK  is  a  polyaromatic  resin,  its  continuous  fiber 
reinforced  composites  have  become  known  as  aromatic  polymer 


composites  (APC). 


Several  different  teams  of  researchers  have  examined 
PEEK  and  APC  properties.  Bishop  compared  the  compressive 
strength  after  impact  behavior  of  an  earlier  version  of  a 
PEEK  composite  (APC1)  with  that  of  carbon  fiber/epoxy  lami¬ 
nates.  She  found  APC1  had  a  higher  compressive  strength 
after  impact  although  exhibiting  a  lower  strength  before 
(1:295) . 

Stober,  et  al ,  confirmed  the  general  solvent  resistance 
of  APC1 ,  and  performed  dynamic  mechanical  analysis  (d.m.a.), 
similar  to  that  performed  in  this  research,  on  both  neat 
resin  and  single  ply  reinforced  composites  (2:1845).  They 
found  "a  strong  correlation"  between  the  a  glass  transition 
and  the  melting  transition  for  neat  PEEK  in  tension  and  APC1 
d.m.a.  in  flexure  (2:1851). 

Jones,  et  al,  outlined  the  mechanical  properties  fur¬ 
ther,  describing  stiffness,  toughness,  and  strength 
(3:1385).  They  confirmed  that  PEEK  composites  deformed  in  a 
nonlinear  fashion  (as  do  thermoplastics  as  a  class),  but 
argued  this  behavior  followed  "a  simple  pattern"  and  that 
small  strain  levels  (  <1%  )  produced  predictable,  linear 
behavior  (3:1386).  This  behavior  has  been  reported  in 
glass-fiber/epoxy  matrix  composites,  too  (21:59)  (23:378). 

A  more  recent  study  by  Cebe,  et  al,  examined  the  ten¬ 
sile  properties  of  PEEK  films  prepared  with  different  ther¬ 
mal  histories  (4:487).  Yield  strains  for  air-cooled  speci¬ 
mens  (processed  at  a  low  quench  rate  similar  to  that  of  the 


neat  resin  sample  employed  in  this  research)  were  5.5%  and 
9.0%  for  test  temperatures  of  125°C  and  25°C,  respectively 
(4:495) . 

Tung  and  Dynes  characterized  the  effect  of  processing 
conditions  on  the  morphology  of  APC2,  the  follow-on  to  APC1 
(5:505).  They  employed  torsion  of  rectangular  samples  in 
their  d.m.a.  over  the  temperature  range  from  -120° C  to 
350 °C,  but  did  not  vary  the  dynamic  frequency  (5:512).  The 
analysis  of  unidirect ional  lamina  in  flexure  has  some  advan¬ 
tages  over  torsion  of  rectangular  samples.  A  more  unified 
procedure  available  for  characterizing  the  mechanical  prop¬ 
erties  of  fiber  reinforced  composites  has  been  shown  to  be 
the  torsion  of  thin  walled  tubes  (27:37).  However,  the 
rectangular  specimen  is  easier  to  fabricate,  and  flexure 
test  is  less  difficult  to  control  (28:801).  Thus  the 
flexure  test  is  more  convenient. 

This  study  attempts  to  predict  the  complex  flexural 
modulus,  E  (u»),  for  an  anisotropic  laminated  beam,  as  a 
special  case  of  laminated  plate  theory.  Experiments  deter¬ 
mined  the  viscoelastic  behavior  of  a  PEEK/graphite  fiber 
composite,  available  commercially  as  APC-2.  Dynamic  three 
point  flexure  tests  measured  the  stiffness  of  thin,  unidi¬ 
rectional  beams.  The  results  from  the  analysis  are  compared 
to  experimental  data.  Also,  f requency- temperature  shifting 
factors  obtained  from  experiments  are  compared  to  values 
predicted  from  the  Williams-Landel-Ferry  (WLF)  equation. 


The  Dynastat  viscoelastic  static  and  dynamic  analyzer 


was  used  for  the  flexural  tests.  It  calculates  the  loads, 
complex  viscoelastic  moduli,  and,  through  a  curve-fitting 
algorithm,  obtains  the  frequency-temperature  shift  factors. 
It  has  a  data  processing  system  that  produced  many  of  the 
graphs  in  this  report.  The  Dynastat  was  designed  to  trans¬ 
versely  flex  long,  thin  samples  of  extremely  stiff  mate¬ 
rials.  This  stiffness  might  otherwise  complicate  tradi¬ 
tional  tensile  experiments  due  to  the  tensile  machine  com¬ 
pliance  approaching  the  material's  compliance.  Typical 
flexure  specimens  in  this  study  are  rectangular  beams  ap¬ 
proximately  42mm  long,  12mm  wide,  and  1mm  thick. 

Although  the  Dynastat  measures  dynamic  flexural  visco¬ 
elastic  properties  with  some  precision,  caution  must  still 
be  used  when  applying  flexural  moduli  to  design  problems. 

The  effects  of  shearing  stress  in  composite  laminated  struc¬ 
tures  can  be  significant  and  must  therefore  be  examined 
before  a  comparison  of  flexural  and  tensile  moduli  can  be 
made.  However,  when  the  structure  is  intended  to  sustain 
bending  moments,  properties  obtained  from  flex  tests  are  of 
value.  The  design  of  structures  from  polymeric  composites 
which  are  intended  to  dampen  vibrations  also  requires  vis¬ 
coelastic  moduli. 

The  goal  of  this  research  was  to  study  the  effect  of 
matrix  viscoelasticity  on  the  mechanical  behavior  of  the 
PEEK  composites,  extending  the  previous  research  just  pre- 
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sented.  Conversely,  what  effects  did  the  presence  of  the 
fibers  have  on  the  viscoelastic  behavior?  A  dynamic 
flexure,  rather  than  tensile,  test  was  selected  to  achieve 
this  goal  because  of  the  high  sample  rigidity,  the  need  for 
high  machine  sensitivity,  and  the  consequent  reduction  in 
the  applied  load. 

One  of  three  objectives  of  this  research  is  to  evaluate 
the  flexural  test  and  to  account  for  other  effects,  i.e., 
shear  and  Poisson  effects.  The  chosen  analytic  model  ap¬ 
plies  a  hybrid  micro-  and  macromechanical  approach  to  the 
analysis  of  a  flexure  specimen.  An  isotropic  sample  of  PEEK 
resin  provided  baseline  data  for  the  matrix  properties.  By 
substituting  that  frequency-dependent  viscoelastic  data  into 
relations  ordinarily  used  for  predicting  the  static  behavior 
of  homogeneous,  anisotropic  elastic  materials,  it  was  hoped 
to  obtain  values  that  agreed  with  experimental  results. 

This  model  is  developed  in  the  second  chapter  of  this  re¬ 
port  . 

The  second  objective  of  this  study  was  to  characterize 
the  frequency  and  temperature  dependency  of  APC-2  in  dynamic 
flexure.  Factors  considered  were  the  longitudinal  and 
transverse  orientations  of  the  continuous  fibers  relative  to 
the  longitudinal  axis  of  the  flexure  specimen. 

Processing  differences,  such  as  quench  rate  during 
specimen  preparation,  were  not  variables.  A  single  unidi¬ 
rectional  plate  provided  the  composite  specimens  for  experi- 


mentation.  The  mid-span  displacements  were  small,  intended 
to  keep  curvatures  small  and  the  strains  in  a  linear  range. 
The  third  section  o£  this  report  details  the  experimental 
procedures  used  to  determine  the  material  behavior. 

A  final  objective  was  to  compare  frequency-temperature 
shift  factors  obtained  in  three  ways:  analytically,  using 
the  WLF  equation;  experimentally,  using  the  Dynastat  curve¬ 
fitting  algorithm;  and  manually,  using  hand-shifted  values 
from  experimental  results.  This  comparison,  and  other  re¬ 
sults  of  the  experimental  work,  are  discussed  in  the  fourth 
chapter  of  this  report. 


I I .  Theory 

Investigations  of  the  dynamic  mechanical  properties  of 
polymers  have  employed  torsion,  tension,  and  three-point 
bending.  Flexure  of  homogeneous  viscoelastic  specimens  has 
provided  reliable  data  that  correlated  well  with  tensile 
properties  (24:9).  Because  the  sample  stiffness  may  ap¬ 
proach  that  of  the  structural  parts  of  the  testing  machine, 
the  experimenter  who  determines  the  viscoelastic  behavior  of 
a  stiff  composite  will  often  choose  flexure  as  the  test 
method.  The  analysis  of  the  experiment  is  then  complicated 
by  the  specimen  geometry,  loading  and  support  conditions,  as 
well  as  any  material  anisotropy  or  nonlinearity.  For  these 
reasons,  bending  tests  of  composites  materials  are  typically 
limited  to  quality  control  and  materials  acceptance  pur¬ 
poses  . 

This  research  attempts  to  model  the  quasi-static  behav¬ 
ior  of  a  dynamic  flexure  test  and  account  for  complications 
due  to  the  laminated  structure  of  the  composite  and  to  the 
material  anisotropy.  APC-2  material  nonlinearity  has  al¬ 
ready  been  noted  by  Jones,  but  small  strains  will  be  assumed 
in  this  analysis  and  enforced  in  the  experimentation. 

A  precise  viscoelastic  model  will  not  be  developed, 
i.e.,  time-dependent  constitutive  relations  will  not  be 
used,  nor  will  the  elastic  solution  for  the  beam  bending 
problem  be  transformed  to  the  Laplace  domain.  Inversion 


from  the  Laplace  domain  will  thus  not  be  undertaken  to 


develop  time-dependent  solutions,  because  the  material  be¬ 
havior  in  the  frequency  domain  is  of  interest. 

Instead,  inertial  effects  will  be  ignored  and  fre¬ 
quency-dependent  material  parameters  will  be  substituted 
into  the  quasi-static  elastic  beam  solution.  The  model  used 
in  this  study  is  built  on  the  analysis  of  the  flexure  test 
for  laminated  composite  materials,  proposed  by  Whitney, 
Browning,  and  Mair  (6:30).  Although  their  procedure  matched 
tensile  properties  quite  well  with  flexural  characteristics 
for  unidirectional  laminates,  they  warned  of  its  injudicious 
use  for  general  laminates  (6:42).  Some  factors  they  consid¬ 
ered  were  free-edge  effects,  transverse  shear  deformation, 
and  shear  coupling.  These  factors  will  be  discussed  further 
below.  Their  findings  form  a  part  of  the  basis  of  the 
American  Society  for  Testing  and  Materials  ( ASTM )  standard 
test  method  for  flexural  properties  of  unreinforced  and 
reinforced  plastics  (7:410). 

A  hybrid  micro-  and  macromechanical  approach  suggested 
by  Sims  and  Halpin,  and  dependent  upon  the  work  of  Halpin 
and  Tsai,  has  shown  that  only  one  viscoelastic  parameter 
need  be  measured  to  predict  composite  behavior  in  the  glassy 
region  of  the  resin  (11:55).  However,  they  assumed  a  time- 
independent  Poisson's  ratio  for  the  resin  below  the  glass 
transition  temperature  (Tg  )  of  the  resin.  Further  develop¬ 
ment  of  their  approach  is  required  for  the  flexural  response 
of  laminated  composites  above  Tg  . 


This  research  will  pre- 


diet  APC-2  flexural  properties  based  on  matrix  and  fiber 
frequency-dependent  behavior  in  a  manner  similar  to  the  Sims 
and  Halpin  rule  of  mixtures. 

Figure  1  depicts  the  flexure  and  laminate  geometry  used 
in  this  research.  According  to  Whitney,  et  al,  Hoff  showed 

that  the  beam's  equivalent  stiffness,  E-.  ,  I,  could  be  attri- 

b 

buted  to  the  bending  stiffness  of  each  lamina,  summed  over 
the  thickness  (6:31): 


where 


n 


■Z« 


(k)  (k) 

11  1 


(1) 


E. ,  =  effective  bending  modulus 

1 1  b 

I  =  moment  of  inertia  of  the  beam  relative  to  the 
midplane 

n  =  number  of  laminate  layers 

(k) 

E^  =  modulus  of  kth.  layer  with  respect  to  the  beam 

axis 

(k) 

I  =  moment  of  inertia  of  kUi  layer  with  respect  to 
the  midplane 

This  relation  holds  for  layered  beams  in  which  the 
plies  are  oriented  symmetrically  about  the  midplane  and  the 
orthotropic  axes  of  material  symmetry  in  each  ply  are  paral¬ 
lel  to  the  beam  edges.  Specially  orthotropic  unidirectional 
beams  are  then  treated  as  macroscopically  homogeneous. 


Whitney,  Browning,  and  Mair  continued  the  approach  that 
Hoff  took,  and  developed  the  constitutive  relations  in  bend- 


ing.  Starting  with  the  constitutive  relation  for  a  lami¬ 
nated  anisotropic  plate  in  pure  bending,  assuming  no  in¬ 
plane  stress  resultants: 
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where  M  ,  M  ,  and  M  are  the  resultant  bending  and  twisting 

x  y  xy 

moments  per  unit  length.  The  plate  bending  and  twisting 

curvatures  are  K  ,  K  ,  and  K  .  In  unidirectional  compos- 

x  y  xy 

ites,  where  the  fiber  direction  is  aligned  with  a  plate 
principal  axis,  the  plate  bending  stiffness  matrix,  [ D ] ,  is 
not  fully  populated.  Because  there  are  no  off-axis  plies, 
D16  =  d26  =  °*  The  inversion  of  the  moment -curvature  rela¬ 
tion  is 


[  D  "  ] 


Mx 

My 
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(3) 


The  D  ”  and  D."  elements  of  the  inverse  of  the  stiffness 

1 2  i  o 

matrix,  ID"],  are  associated  in  name  with  the  Poisson  effect 
and  shear  coupling,  respectively. 

To  neglect  the  effects  of  shear  coupling,  one  continues 
by  assuming: 


My  »  Mxy  "  0 


(4  ) 
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and  that  the  vertical  deflection,  w,  is  a  function  of  x, 
i.e.,  w  =  w(x).  This  is  applicable  for  beams  which  have  a 
high  length-to-width  ratio. 

As  Whitney,  et  al,  point  out,  this  approach  neglects 
shear  coupling  The  effect  of  shear  coupling  may  not  neces¬ 
sarily  be  ignored  in  angle-ply  laminates,  and  the  dependence 
of  the  deflection  along  the  y-direction,  or  width,  will 
cause  warping  at  the  beam  supports.  This  research  is  con¬ 
ducted  only  on  unidirectional  beams,  and  their  length-to- 
width  ratio  is  quite  high. 

The  dependence  of  the  deflection  on  the  width  results 
from  the  expressions  for  the  curvatures  of  a  plate  in  terms 
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Ref.  8,  the  relation  which  accounts  for  shear  coupling 
between  the  experimentally  determined  effective  flexural 


V' 


modulus,  E ^  ,  and  the  Young's  modulus  in  the  1-direction, 

b 

E^,  and  which  produced  reasonably  good  experimental  re- 

(L/b)2E  ^ 

<7> 

(L/b)  *•  -  V  {E  /E  ) 

12  11  22 

=  Young's  moduli  in  the  1-  and  2-  directions, 
respectively 
=  span-to-width  ratio 

*  Poisson's  ratio  for  the  strain  in  the 
2-direction  when  the  beam  is  stressed  in 
the  1-direction 
The  A3TM  recommendation  for  the  value  of  L/b  is  at  least 
1.6,  for  the  three-point  loading  of  a  reinforced  plastic 
flexural  specimen  with  a  length-to-thickness  ratio  of  40 
(7:417) . 

Following  the  approach  Ashton  and  Whitney  took  (Refer¬ 
ence  22),  and  that  taken  by  Schultz  and  Tsai  (23:369),  Whit¬ 
ney,  et  al  (6:34),  expressed  the  constitutive  relations  for 
a  homogeneous  orthotropic  beam  in  terms  of  the  stresses  and 
strains  as: 


1 1 

suits,  is: 
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where 


E  ,  E 
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L/b 
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where  the  Qij  are  the  moduli.  This  approach  is  identical  to 
Tsai's  transformation  of  the  full  stress  tensor  because 
there  are  no  off-axis  plies  (25:440).  From  an  equilibrium 
study  of  a  laminated  plate: 


dM 


-  Q  =0 

x 


where  the  total  transverse  shear  resultant,  Q,  is  Q  =  bQ  , 

x 

or,  the  product  of  the  plate  (beam)  width  and  the  transverse 
shear  resultant  per  unit  width.  Because  of  the  assumption 
in  Eqn .  (4),  Q  =  dM/dx,  where  M  is  the  applied  bending 
moment,  and  is  the  product  of  the  resultant  bending  moment 
per  unit  width  and  the  plate  (beam)  width.  It  is  now  pos¬ 
sible  to  account  for  the  rotation  of  beam  elements  due  to 
shear . 

Denoting  the  normals  to  the  midplane  of  the  beam  rela¬ 
tive  to  the  x  and  the  y  axes,  respectively,  as  &  ,  and  w  , 

x  y 

the  constitutive  equations  for  plate  bending  are  analogous 
to  the  Timoshenko  beam  theory.  Changing  Whitney's  sign 
convention  to  that  used  by  Shames  and  Dym  (30:342),  one 
obtains  the  curvature-moment  relation  similar  to  Eqn.  (3): 


:v  o 


w 


»♦  V  V 


result  into  Eqn.  (10)  produces  the  governing  differential 
equation  for  the  elastic  flexure  of  a  beam: 


« 

p 


‘v-‘ 


d2v/dx2  =  ( 1/kGbh ) ( d2  M/dx2 )  -  M/EI 


(13) 


With  the  applied  moment  known  as  a  function  of  x,  the  elas¬ 
tic  deflection  can  be  obtained  from  the  integration  and 
proper  use  of  boundary  conditions  in  Eqn.  (13). 

For  the  s imply-suppor ted  beam  under  midspan  deflection. 


the  maximum  deflection,  w  ,  occurs  at  the  midspan,  and  is 


PL 


PL 


w  = 

o 


(14) 


4kGbh  4Ebh 


where  the  applied  load  at  midspan  is  P.  Rearranging  terms 
for  the  effective  flexural  modulus  yields: 


PL 


E  = 


(15a) 


PLh 


(  4bh  v  - 


Gk 


The  net  effect  of  shear  deformation  is  to  decrease  the 
calculated  value  of  the  bending  modulus  if  shear  were  ig¬ 
nored.  That  part  of  the  deflection  due  to  shear  can  not  be 
measured  directly  in  a  flexure  test.  To  calculate  the 
contribution  of  shear,  and  attempt  to  find  the  part  of  the 
effective  flexural  modulus  due  only  to  bending,  the  measured 
values  of  the  effective  flexural  modulus  will  be  increased 
by  including  the  negative  term  in  the  denominator  in  Eqn. 
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(15a)  to  the  portion  of  the  denominator  without  the  shear 
term.  The  experimental  equipment  in  this  research  does  not 
account  for  shear  in  its  automatic  scaling  of  raw  data,  and 
instead  calculates  the  effective  flexural  modulus  as 

E  =  PL3/4bh3wQ  (15b) 

For  highly  anisotropic  materials,  the  contribution  of  the 
second  term  in  the  denominator  of  Eqn .  (15a)  may  be  consid¬ 
erable.  The  effects  of  this  contribution  are  examined  in 
this  research. 

Additionally,  free  edge  effects  are  considered  in  the 
model  because  the  state  of  stress  is  actually  three  dimen¬ 
sional  in  a  small  zone  near  the  sides  of  the  beam.  The 
stress  field  can't  be  considered  a  plane  state  for  a  dis¬ 
tance  equal  to  the  thickness  in  from  the  edges.  To  minimize 
the  free  edge  effects,  i.e.,  to  approximate  a  state  of  plane 
stress  across  the  majority  of  the  width  of  the  beam,  the 
flexure  specimen  must  have  a  high  width-to-thickness  ratio, 
b/h.  ASTM  standards  recommend  b/h  =  25  for  a  specimen  with 
a  length-to-thickness  ratio  of  L/h  =  40  (7:417),  when  that 
sample  is  tested  to  fail ure . 

The  maximum  dimensions  of  the  Dynastat  test  chamber, 
however,  forced  some  experimental  trade-offs  in  terms  of 
specimen  geometry.  These  are  discussed  in  the  next  chapter. 
It  suffices  to  say,  that  in  order  to  obtain  measurable  re¬ 
sults  in  a  regime  where  small  angle  approximations  and 


linear  theory  are  still  valid,  the  nominal  b/h  value  for  all 
flexure  specimens  in  this  research  is  12.5.  Loading  did  not 
approach  failure  loads,  and  edge  delaminations  were  not 
observed . 

Up  to  this  point  in  the  analysis,  the  microscopic 
nature  of  the  composite  has  been  suppressed.  To  continue, 
one  assumes  the  matrix  of  the  composite  to  act  isotropical¬ 
ly,  while  allowing  the  fibers  to  behave  anisotropically . 

This  approach,  taken  by  Halpin  and  Tsai  and  reported  in  Ref. 
9,  has  been  called  a  rule  of  mixtures  approach.  The  macro¬ 
scopic  composite  behavior  is  assumed  to  be  described  by  the 
volume-weighted  sum  of  the  behaviors  of  the  fiber  and  matrix 
separately.  Referring  to  Figure  1,  then, 

En  =  k(Ef  V  f  +  E  V  )  (16) 

L  i  t  t  mm 

where 

Ejj  =  Young's  modulus  in  the  1-  direction 

E  c,  E  =  Young's  moduli  of  the  fiber  and  matrix, 

f  m 

respectively 

V  ,  V  =  volume  fraction  of  fiber  and  matrix,  respec- 

t  m 

t i ve 1 y 

k  =  fiber  alignment  factor 

The  fiber  alignment  factor  is  assumed  to  be  0.9  throughout 
this  analysis.  Since  the  prepreg  plies  were  cut  and  stacked 
up  by  hand,  and  a  value  of  k  =  1.0  would  describe  perfectly 
aligned  fibers,  the  above  factor  yields  a  conservative 


In  the  transverse  direction,  the  composite  modulus 


becomes 


Ef  <1  +  *  vf  )  ♦  e Em  (i  -  vf  ) 

E22  *  E,  -  (17) 

Ef  (1  -  Vf  )  +i  Em  (!  +  Vf  /<) 

where  f  is  the  fiber  packing  factor.  This  value  depends  on 
the  spacing  around  individual  fibers  and  their  neighboring 
fibers.  Whitney  has  shown  that  the  transverse,  not  longitu¬ 
dinal,  fiber  modulus  should  be  used  for  Ef  in  Eqn.  (17)  if 
significant  fiber  anisotropy  is  present  (17:188).  If  the 
circular  fibers  are  arrayed  in  a  square,  Hewitt  and  deMal- 
herbe  (10:118)  suggested  calculating  $  as 

£  =  1  +  40V  f10  (18) 

The  rule  of  mixtures  applied  to  the  calculation  of  the 
Poisson's  ratio  yields: 
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=  V 


m  m 


+  V 


(19) 


Of  course.  Equations  16  through  19  make  no  allowance 
for  the  frequency  dependency  of  any  of  the  material  parame¬ 
ters.  This  research  will  employ  Eqns .  (16)  and  (17)  in  a 

quasi-static  manner  with  flexural  properties  substituted  for 
the  Young's  moduli  of  the  fiber  and  matrix  as  appropriate. 
The  effect  of  the  viscoelastic  nature  of  the  matrix  can  then 


be  observed  on  the  applicability  of  the  rule  of  mixtures. 


Assuming  frequency- independent ,  elastic,  anisotropic 
graphite  fiber  properties,  and  an  isotropic  matrix,  only  the 
matrix  material  need  be  modeled  by  frequency-dependent 
variables.  From  Eqns.  6  and  15,  with  shear  deformation 
neglected,  two  material  properties  to  be  described  visco- 
elastically  are  Poisson's  ratio  and  the  flexural  modulus. 
Temperature  dependency  will  observed  and  related  to  fre¬ 
quency  by  superposition. 

The  frequency  dependent  Poisson's  ratio  and  flexural 
modulus  will  consist  of  in-phase  and  out-of-phase  compo¬ 
nents.  However,  an  extensive  study  by  Koppelmann,  reported 
by  Ferry,  indicates  that  "the  imaginary  component  of  Pois¬ 
son's  ratio  ...  is  small,  perhaps  zero"  (11:444).  Then  the 
matrix  Poisson's  ratio  is  ym  -  ym  (o»)  ~Vm  '  iu)  • 

For  the  harmonic  strain  excitation  to  be  considered 
here,  the  complex  modulus,  E  (cj)  ,  is  defined  as  the  ratio  of 
the  harmonic  stress  response  to  the  strain.  If  sufficient 
time  after  the  application  of  strain  has  occurred  to  neglect 

the  initial  conditions,  for  strain  €(t)  =e  explicit],  and 

o 

stress  a  ( t )  =  a  explicit  1, 

o 

E  =  a  /€  =  E  *(<J)  =  E'(C*»)  +  iE”(w)  (20) 

o  o 

E'(w)  is  the  frequency  dependent  in-phase  component  of  E  (cj) 
(20:  24  ).  E'(£d)  is  also  called  either  the  storage  modulus  or 

the  real  part  of  the  complex  modulus.  E'’(Ct>)  is  called 
either  the  out-of-phase,  or  imaginary,  component,  or  the 


loss  modulus.  The  ratio  of  the  out-of -phase  to  in-phase 
components  is  called  the  loss  factor.  The  loss  factor  is 
the  tangent  of  the  phase  angle, 6  ,  between  the  components  of 
the  complex  modulus,  i.e.,  E'VE'  =  tan  d  .  The  loss  factor 
is  related  to  the  material  damping  factor  (21:56). 

When  the  fibers  of  a  unidirectional  composite  beam  are 
oriented  such  that  the  fibers  are  parallel  with  the  longitu¬ 
dinal  axis,  dynamic  three  point  bending  determines  E ,  ^ "  (ref¬ 
erence  Figure  1).  Similarly,  when  the  fibers  are  oriented 
such  that  they  are  transverse  to  the  beam  longitudinal  axis, 
dynamic  three-point  bending  determines  E^". 

The  important  characteristic  of  materials  modeled  as 
thermorheologically  simple  solids  is  their  frequency-temper¬ 
ature  behavior.  The  flexural  modulus  at  one  temperature 
will  exhibit  a  frequency  dependence  that  graphically  matches 
a  portion  of  the  behavior  at  another  temperature.  Finite 
temperature  changes  produce  shifts  in  the  modulus  on  a  loga¬ 
rithmic  frequency  scale.  The  Williams,  Landel,  and  Ferry 
(WLF)  equation  predicts  the  amount  of  this  shift  as  a  func¬ 
tion  of  temperature.  Experimentally,  one  obtains  the  master 
curve  of  a  material  by  shifting  spectra  over  the  actual 
frequency  range  relative  to  a  reference  temperature.  Thus 
the  behavior  of  the  modulus  can  be  presented  relative  to  a 
reduced  frequency  over  a  wider  range  of  frequencies  than 


The  reduced  frequency  is  the  product  of  the  actual 
frequency  (at  which  modulus  data  was  obtained)  and  the 
amount  of  shift,  or  shift  factor.  A  good  approximation  of 
the  shift  factor,  a(T),  for  a  reference  temperature,  T 

K 

(about  50°C  above  the  glass  transition)  that  holds  for  many 
polymers  is 

log  a (T)  =  - 8 . 86 ( T  -  T_)/(101.6  +  T  -  T  )  (21) 

1 U  K  K 

where  T  is  the  actual  temperature  of  the  measured  spectra. 
Eqn.  (21)  is  the  WLF  equation.  It  holds  only  near  and  above 
the  glass  transition  temperature,  T  (11:289). 

Quite  apart  from  the  empirical  curve-fitting  technique 
of  horizontally  shifting  the  spectra  (or  the  application  of 
the  WLF  equation)  is  another  temperature  reduction  on  the 
complex  modulus.  This  is  predicted  by  Rouse's  theory.  Fer¬ 
ry  reports  that  this  is  based  on  the  hypothesis  "that  the 
modulus  contributions  are  proportional  topT"  (11:269).  The 
density,  p  ,  is  often  dropped.  This  is  acceptable  because 
"it  is  impossible  to  distinguish  empirically  between  the 
merits  of  pT  and  T"  (11:271). 

Master  curves  generated  in  the  results  section  of  this 
report  plot  the  moduli,  multiplied  by  the  factor  of  T  /T, 
versus  the  reduced  frequency,  f  x  a(T),  where  f  is  in  hertz. 


III.  Experimental  Procedure 

Mater ia Is 

Specimens  in  this  research  were  cut  from  a  6  x  9  inch 
composite  plate  prepared  by  hand  from  multiple  plies  of  PEEK 
tape  pre-impregnated  with  graphite  fiber.  This  tape,  or 
prepreg,  was  manufactured  by  the  Fiberite  Co.,  a  subsidiary 
of  I  Cl  Americas,  Wilmington,  Delaware,  and  was  designated 
APC-2/AS4 .  Individual  plies  were  0.005mm  thick.  Eight  ply 
laminates  were  prepared  in  a  unidirectional  stacking 
sequence,  and  processed  under  temperature  and  pressure  to 
produce  composite  specimens  from  1.016mm  to  1.040mm  thick. 
The  processing  procedure  described  below  results  in  compos¬ 
ites  with  properties  similar  to  those  listed  in  Table  1. 

Table  1 

Typical  Mechanical  Properties  of  APC-2  at  25°C 
Data  from  (12:1410)  and  (3:1388) 


Flexural  Modulus,  GPa 

134 

Longitudinal  Tensile  Modulus,  GPa 

124 

-  141 

Transverse  Tensile  Modulus,  GPa 

8.9  - 

10 . 3 

In-Plane  Shear  Modulus,  GPa 

4.0  - 

6.27 

V 

0.38 

L  2 

V21 

** 

o 

o 

Melting  Point,  °C 

345 

The  composite  processing  followed  the  thermal  history 
optimized  by  the  Fiber ite  Co.  (13),  and  was  conducted  at  the 
Air  Force  Materials  Laboratory  Non-Metall ics  Division, 
Structural  Materials  Branch.  The  prepreg  plies  were  placed 
in  a  mold  press  under  a  light,  contact  pressure.  The  press 
platens  were  heated  to  385° C  from  room  temperature  in  60 
minutes.  Platen  pressure  was  then  increased  to  200  psi,  and 
maintained  for  45  minutes  at  388. 1° C  +  /-  5.0°C.  A  slow 
quench  rate  period  followed,  still  under  high  pressure,  from 
391°C  to  46°C  for  85  minutes. 

The  average  quench  rate  was  -4°C/min.  This  produces  a 
degree  of  crystall inity,  by  density,  of  about  .45  -  .48 
(5:511).  The  degree  of  crystallinity  is  a  measure  of  order 
of  the  polymeric  matrix,  with  an  amorphous  state  having  zero 
crystallinity.  Nonlinear  viscoelastic  behavior  has  been 
observed  with  increasing  crystallinity.  Other  measured 
properties  are  listed  in  Table  2.  An  acid  digestion  tech¬ 
nique  determined  the  volume  fractions. 

Table  2 

APC-2  Morphology 


Specific  Gravity  1.607 

Fiber  Volume  Fraction  0.6401 


Matrix  Volume  Fraction 


0 . 3465 


Quality  control  measures  included  a  microscopic  study 
of  a  sample  cross-sectional  area  of  the  plate,  at  magnifica¬ 
tions  up  to  lOOOx,  and  an  ultrasonic  scan  at  3.5dB.  No 
significant  internal  voids  or  delaminations  were  observed, 
although  the  sum  of  the  volume  fractions  indicates  a  void 
content  of  less  than  2%. 

Flexure  samples  were  cut  from  the  plate  to  nominal 
dimensions  of  2  x  1/2  inches.  Samples  with  the  fibers 
aligned  with  the  length  are  termed  "0°"  specimens,  while 
"90°"  beams  have  the  fibers  oriented  with  the  width.  A 
water-cooled,  diamond  blade  saw  cut  all  specimens,  and  some 
edges  were  wet  sanded  to  final  shape,  but  only  along  the 
lengthwise  direction,  in  accordance  with  ASTM  D  790M 
(7:411).  Neither  the  upper  nor  lower  surfaces  were  ma¬ 
chined.  The  orientations  of  the  0°  and  90°  samples  with 
respect  to  the  three  point  flex  test  are  shown  in  Figure  1. 

Neat  resin  samples  were  prepared  from  molding  several 
layers  of  5  mil  thick  PEEK  film  in  a  high  temperature  press. 
Temperatures  and  quench  rates  were  nearly  identical  to  those 
for  the  composite  laminate,  although  the  pressure  was  only 
50  psi  at  the  high  temperature.  A  nominal  specimen  thick¬ 
ness  of  2mm  was  obtained.  A  typical  Young's  modulus  value 
of  2.6  GPa  was  supplied  by  the  manufacturer,  ICI,  for  the 
non-crystallized  Stabar  K  film  (14:1).  Cebe,  et  al,  have 
reported  mechanical  properties  of  crystallized  PEEK  films 
from  0.005cm  to  0.010cm  thick.  For  a  similar  thermal  his- 


tory  as  in  this  research,  their  values  for  Young's  modulus 
were  from  1.4  GPa  to  2.0  GPa,  with  a  degree  of  crystallinity 
of  about  0.32  (4:490).  Differences  between  composite  and 
neat  resin  behavior  were  thus  expected  to  occur  due  to  the 
contributions  of  the  fiber  to  the  rigidity  and  to  crystal 
formation . 

Both  the  Dynastat  flexure  supports  and  its  environmen¬ 
tal  chamber  determine  the  maximum  allowable  specimen  widths 
and  lengths.  For  three-point  bending,  the  maximum  specimen 
width  is  0.5  inches.  The  inner  diameter  of  the  environmen¬ 
tal  chamber  is  approximately  3  inches.  The  span  of  the 
lower  support  is  41.71  mm,  so  most  specimens  were  sized  just 
above  this  length. 

For  highly  anisotropic  laminates,  the  minimum  length- 
to-thickness  ratio  should  be  at  least  40,  with  60  recom¬ 
mended  by  the  ASTM  in  reference  7.  Shear  deformation  other¬ 
wise  will  be  significant.  For  the  Dynastat  support  span, 
this  higher  value  prescribes  a  sample  thickness  below  1mm. 

A  standard  sample  geometry  is  desirable  for  both  fiber 
orientations  and  for  the  neat  resin  flexure  specimen.  It  is 
felt  that  a  sample  thickness  less  than  1mm  would  lower  the 
stiffness  of  the  90°  and  neat  resin  specimens  to  the  point 
that  the  load  cell  would  have  difficulty  distinguishing 
force  levels  from  noise  during  high  temperature  flexure. 
Another  desirable  test  parameter  is  a  small  maximum  center 
span  deflection  that  produces  a  small  maximum  outer  fiber 
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strain.  This  would  not  be  possible  with  extremely  thin, 
flexible  specimens.  Thus  a  trade-off  was  made.  A  rela¬ 
tively  low  width-to-thickness  ratio,  12.5,  was  used,  but  the 
length-to-thickness  ratio  of  40  was  thought  to  be  accept¬ 
able.  The  length-to-width  ratio  of  3.33  is  double  the 
recommended  value  of  1.6  and  further  minimizes  the  effect  of 
shear  coupling.  This  effect  is  calculated  and  presented  in 
the  results  chapter. 

Prior  to  testing,  all  specimens  were  stored  in  dehumi¬ 
difiers.  The  size  of  each  specimen  was  measured  just  prior 
to  testing,  in  accordance  with  ASTM  790M-84  specimen  condi¬ 
tioning  standard  8.1.  These  geometric  values  were  necessary 
inputs  to  Dynastat  automatic  scaling  functions  described 
below. 


Test  Equipment 

The  viscoelastic  character ization  studies  were  per¬ 
formed  on  an  instrument  commercially  available  from  Imass, 
Inc.,  Hingham,  Massachusetts,  and  known  as  the  Dynastat. 
The  system  is  fully  automated  for  data  acquisition  and 
reduction.  The  five  main  components  of  the  Dynastat  (the 
servo  controller,  the  Dynatherm,  the  Dynalyzer,  the  Data- 
lyzer,  and  the  test  frame)  are  illustrated  functionally  in 
Figure  2,  and  described  below. 


1 


The  test  frame  consists  of  an  environmental  chamber,  a 
vertical  push-pull  rod  with  specimen  supports,  two  displace- 
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merit  transducers,  and  a  load  cell.  In  this  study,  the  low 
range  displacement  transducer  provides  both  measurement  of 
and  control  of  the  vertical  displacement  of  the  lower  sup¬ 
port.  The  transducers  are  below  the  sample,  outside  the 
chamber.  A  fixed  indentor  is  positioned  above  the  flexure 
sample.  The  load  cell  is  above  the  sample,  and  is  outside 
the  chamber . 

All  flexure  tests  were  performed  using  three  point 
loading.  The  static  indentor  and  the  moving  lower  supports 
have  a  load  point  diameter  of  5.0  mm. 

The  push-pull  rod  is  rigidly  coupled  to  the  coil  of  a 
linear  motor  underneath  the  chamber,  and  is  mounted  in  an 
air  bearing  for  transverse  stiffness.  The  linear  motor 
consists  of  a  permanent  magnet  and  a  spider  suspended  coil 
wound  on  a  lightweight  coil  form  (15:133). 

The  controlling  servos  for  the  linear  motor  are  located 
in  the  sub-unit  known  as  the  Dynastat  (the  manufacturer 
refers  to  the  entire  test  ensemble  by  that  same  name.  Un¬ 
less  specific  mention  is  made  to  the  servo  controller, 
further  references  to  the  Dynastat  include  the  ensemble). 

The  servo  system  provides  closed-loop  control  using  the 
displacement  transducer  for  feedback  (15:135).  The  push¬ 
button  potentiometers  on  the  servo  controller  "operate  logic 
switches  which  permit  ...  step  function  inputs  of  ...  dis¬ 
placement  and  ...  dynamic  external  inputs"  (16:6).  The 
experimenter  programs  the  static  and  dynamic  amplitudes  into 


the  servo  controller  and  then  includes  timer  inputs,  dynamic 
signals,  and  thermal  triggers  into  the  test  by  adjusting  the 
other  sub-units  of  the  ensemble. 


Temperature  control  of  the  environmental  chamber  is 
provided  by  the  Dynatherm.  The  Dynatherm  also  programs 
temperatures  increments  and  provides  thermal  triggers  to  the 
servo  controller.  Thumbwheel  settings  on  the  front  of  this 
sub-unit  define  the  lower  and  upper  temperature  setpoints, 
and  the  number  of  intervals  between  them.  Up  to  99  inter¬ 
vals  can  be  programmed  into  the  temperature  range  of  between 
-150 °C  and  350°C  (18:2)  . 

The  Dynatherm  controls  the  heating  of  forced  air 
through  the  environmental  chamber  by  sensing  chamber  temper¬ 
ature  with  a  platinum-resistance  thermometer  ( PRT )  and  servo 
controlling  the  heater.  The  control  logic  is  a  mixture  of 
proportional  and  integral  control,  and  provides  temperature 
stability  of  +/-0.1°C.  Typical  response  for  a  programmed 
5° C  increment  is  shown  in  Figure  3.  Figure  3  is  a  portion 
of  a  stripchart  history  recorded  for  a  calibration  run.  The 
vertical  axis  represents  temperature  deviation  from  the 
setpoint.  The  horizontal  axis  represents  time,  illustrating 
that  dynamic  data  collection  sweeps  were  conducted  at  the 
end  of  a  stable  time  soak. 

The  PRT  is  placed  near,  but  not  in  contact  with,  the 
flexure  sample  in  the  chamber.  No  thermocouples  are  at¬ 
tached  to  the  specimen.  The  soak  time  is  intended  to  insure 
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that  the  sample  temperature  is  as  close  as  possible  to  the 
temperature  of  the  chamber  air. 

Dynamic  control  and  signal  correlation  is  provided  by 
the  Dynalyzer.  The  Dynalyzer  and  the  Dynatherm  coopera¬ 
tively  provide  automatic  temperature  and  frequency  sweeps. 
The  Dynalyzer  digitally  generates  a  harmonic  signal  that 
drives  the  linear  motor,  and  which  serves  as  the  reference 
for  in-phase  and  out-of-phase  component  matching.  The  Dyna¬ 
lyzer  obtains  complex  moduli  values  using  a  modified  Fourier 
transform  and  auto-correlation  process  that  compares  the 
load  cell  and  displacement  transducer  outputs  with  the  ref¬ 
erence  signal.  The  frequency  range  of  the  Dynalyzer  is  from 
0.1  Hz  to  199.9  Hz,  with  a  maximum  error  of  +/-  0.016  Hz  in 
the  range  from  0.16  Hz  to  15.9  Hz  used  in  this  research. 

The  range  of  measurement  for  the  loss  factor,  or  tan  6  ,  is 
from  0.001  to  19.999  (19:4). 

The  Dynalyzer  furnishes  four  digital  values  (raw  data 
in  units  of  volts/volts)  to  the  Datalyzer  for  storage  in  a 
data  array.  These  are  the  unsealed  loads  and  displacements, 
and  the  unsealed  in-  and  out-of-phase  components  of  the  load 
relative  to  the  displacement.  The  Datalyzer  is  a  micropro¬ 
cessor-based  instrument  which  contains  a  random  access  mem¬ 
ory  (RAM).  The  read-only  memories  (ROMS)  of  the  Datalyzer 
permit  it  to  communicate  with  peripheral  devices  including  a 
printer,  a  plotter,  and  a  floppy-disc  drive.  The  utility  of 


these  devices  is  discussed  later. 


The  Datalyzer  manipulates  the  raw  data  with  geometry 
factors,  load  and  displacement  scales,  and  several  correc¬ 
tion  factors  such  as  machine  compliance.  These  factors  are 
entered  by  keyboard.  In  displacement  control  mode,  the 
Datalyzer  calculates  E*((j)  from  the  relation  for  an  isotrop¬ 
ic,  homogeneous  beam  without  shear,  i.e.,  from  Eqn .  (15b). 

The  Datalyzer  output  can  be  printed  in  real  time,  or  stored 
on  floppy  diskette  for  later  analysis.  The  ROMS  provide  the 
capability  to  plot  the  results  of  the  dynamic  tests  against 
appropriate  frequency  or  temperature  scales. 

The  Datalyzer  assigns  an  eighth-order  polynomial  to  the 
moduli  data  at  each  temperature  as  a  function  of  the  actual 
frequency.  The  Datalyzer  then  computes  the  shift  factors 
for  each  spectrum  with  respect  to  the  spectra  at  the  next 
higher  and  lower  temperatures.  First  the  Rouse  theory,  or 
temperature  correction,  is  applied  to  the  moduli.  Then  the 
minimum  integral  area  under  adjacent  polynomials  is  calcu¬ 
lated.  Each  calculation  results  in  a(T)  for  that  tempera¬ 
ture  . 

A  real  time  check  on  the  operations  of  the  Dynastat  is 
provided  by  an  oscilloscope.  The  oscilloscope  displays  the 
output  waveforms  of  the  load  and  displacement  channels 
during  frequency  sweeps.  This  function  is  necessary  to 
ensure  the  indentor  remains  on  the  specimen  at  all  times. 
Clipped  sine  waves  on  the  load  trace  indicate  that  the 
indentor  has  i i sen  off  the  specimen,  and  that  the  indentor 


requires  readjustment.  Also,  the  noise  in  either  channel  is 
visually  evident  at  any  time  during  test  runs.  This  visual 
monitoring  of  the  waveforms  was  essential  to  determining  the 
magnitudes  of  the  static  and  dynamic  programs  on  the  servo 
controller.  Later,  the  monitoring  was  essential  to  deter¬ 
mining  the  quality  of  the  resulting  data. 


Experimental  Method 

The  viscoelastic  data  in  this  study  resulted  from  dy¬ 
namic  flexure  only.  However,  before  a  dynamic  test  was 
conducted  on  a  given  ply  orientation  or  on  the  neat  resin 
specimen,  the  linearity  of  the  static  load-displacement 
relation  had  to  be  determined  at  several  temperatures  within 
the  expected  range  of  the  dynamic  test.  In  addition,  an 
optimum  ratio  of  the  static  to  the  dynamic  deflections  had 
to  be  determined.  Obviously,  the  amplitude  of  the  maximum 
dynamic  deflection  had  to  be  less  than  the  static  load,  or 
the  indentor  would  rise  off  the  sample.  But  the  need  for  a 
high  signal  to  noise  ratio  was  the  primary  reason  for  deter¬ 
mining  the  static-to-dynamic  displacement  ratio. 

During  static  load-displacement  tests,  the  sample  was 
loaded  in  flexure,  and  the  displacement  was  stepwise  in¬ 
creased  to  0.1%  strain,  as  calculated  from  the  relation  for 
the  ASTM  standard  : 


e  =  6hw  /L' 
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No  strain  gauges  were  placed  on  the  sample.  Experimen¬ 
tal  data  from  Whitney,  et  al  (6:38),  and  Srinivasan,  et  al 
(26:279),  revealed  little  difference  in  measured  bending 
stiffness  compared  to  the  method  of  displacement  measurement 
used  here.  The  maximum  outer  fiber  strain  of  0.1%  is  well 
below  the  ASTM  recommended  limit  of  5%  due  to  flexure 
(7:411). 

These  incremental  deflections  were  repeated  at  elevated 
temperatures,  and  the  output  of  the  load  cell  and  the  dis¬ 
placement  transducer  was  recorded.  The  results  are  dis¬ 
cussed  in  the  next  chapter. 

The  maximum  dynamic  displacement  was  determined  by  a 
preliminary  dynamic  test,  the  purpose  of  which  was  to  find 
an  optimum  s  ignal -to-no ise  ratio  (S/N)  in  the  load  channel. 
As  the  sample  softened  with  increasing  temperature,  S/N 
could  be  expected  to  worsen,  degrading  the  accuracy  of  the 
auto-correlation  process  of  the  Dynalyzer.  All  dynamic 
mechanical  tests  were  then  conducted  at  this  deflection 
ratio.  This  ratio,  and  its  relation  to  the  static  load- 
displacement  response,  are  discussed  in  the  next  chapter . 

The  temperature/dynamic  mechanical  tests  were  conducted 
after  machine  calibration  and  sample  geometry  factors  were 
entered  into  the  Datalyzer.  This  entry  was  required  only  at 
the  start  of  the  test.  With  the  exception  of  the  neat  resin 
sample,  all  frequency  sweeps  were  conducted  at  discrete 
frequencies  in  a  two  decade  band  from  0.16  Hz  to  15.90  Hz  (1 
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rad/sec  to  100  rad/sec).  Neat  resin  measurements  were  at¬ 
tempted  at  the  same  frequencies,  but  only  the  data  at  four 
discrete  frequencies  between  2.5  Hz  and  10.09  Hz  were  of  any 
value  due  to  noise.  The  sweep  (and  data  collection  period) 
for  all  samples  started  7.0  minutes  after  the  environmental 
chamber  temperature  reached  +  /-  1.0°C  of  the  programmed 
temperature . 

The  chamber  temperature  ramped  up  after  each  frequency 
sweep  in  10°C  steps  for  the  composite  samples,  but  to  only 
three  setpoints  for  the  neat  resin.  No  cryogenic  tests  were 
run,  but  rather  the  temperature  ranges  were  typically  from 
near  ambient  to  about  250°C. 

The  Dynalyzer  averaged  the  Dynastat  output  over  10  to 
15  cycles  at  each  frequency  to  further  increase  noise  dis¬ 
crimination.  An  improper  assumption  that  data  at  the  low 
temperatures,  or  glassy  region,  would  not  be  as  noisy  led  to 
the  use  of  10  cycles  versus  15  cycles  in  that  region.  This 
assumption  resulted  in  noisy  data  in  this  region,  although 
much  of  the  scatter  may  be  attributable  to  composite  behav¬ 
ior. 

After  transmitting  the  averaged  data  to  the  Datalyzer, 
the  Dynalyzer  advanced  to  the  next  driving  frequency,  while 
the  Dynatherm  ensured  chamber  temperature  stability.  The 
Dynalyzer  delayed  accepting  digital  data  from  the  Dynastat 
at  any  new  frequency  for  an  adjustable  time  period.  This 
delay  equaled  one  cycle  for  frequencies  less  than  1.0  Hz, 


and  equaled  1  second  for  frequencies  at  and  above  1.0  Hz. 

The  delay  ensured  transients  had  damped  out  after  the  in¬ 
crease  in  driving  frequency. 

An  entire  frequency  sweep,  conducted  over  fifteen 
cycles  for  eleven  discrete  frequencies  at  one  temperature, 
took  approximately  five  minutes.  Heating  to  the  next  tem¬ 
perature  took  approximately  twelve  minuses  (this  time  in¬ 
creased  at  higher  temperatures).  After  a  seven  minute  soak, 
the  sweep  was  conducted.  Usually  an  eight  hour  period  was 
required  to  ramp  up  to  all  temperatures  and  conduct  sweeps 
at  each  one. 

Post  test  procedures  included  data  transmission  from 
RAM  to  disk  and  scaling  of  the  plotter  axes  from  the  range 
of  recorded  data.  Shift  factors  for  frequency-temperature 


superposition  were  calculated.  The  resultant  master  curves 
were  plotted  for  the  composite  samples  only. 


I V .  Results 


Before  the  major  findings  of  this  research  are  ex¬ 
plored,  the  results  of  the  static  load-d isplacement  tests 
will  be  presented.  Those  results  defined  the  region  of  the 
dynamic  sweep,  that  is,  the  static-to-dynamic  displacement 
ratio.  The  temperature  dependence  of  E'  (cj)  ,  E' '(<*/),  and  the 
loss  factor,  tan  6  ,  will  then  be  presented.  Frequency- 
temperature  shift  factors  are  next  presented  and  compared. 
The  f requency- temperature  transformation  is  only  applicable 
near  or  above  T^  ,  therefore  the  shift  factors  in  this  re¬ 
search  are  examined  only  in  that  temperature  region.  Master 
curves,  generated  by  the  f requency- temperature  superposition 
of  individual  spectra,  are  presented,  and  compared  by  shift 
technique.  Experimental  results  are  compared  to  the  theo¬ 
retical  predictions  for  shear  deformation  effects,  shear 
coupling,  and  volume  fraction  contributions. 

The  results  of  the  static  deflection  tests  are  depicted 
in  Figures  4  and  5.  Figure  4  plots  the  load  cell  output  (?) 
as  a  function  of  the  maximum  displacement  (w{))  of  the  push- 
pull  rod  for  the  90°  specimen.  This  composite  demonstrates 
a  linear  P  -  w()  relation  in  a  temperature  range  of 
24.4  °C  <  T t)  <  250.  3  °C.  The  maximum  deflection  correspondr 
roughly  to  a  maximum  outer  lamina  strain  of  0.1%,  as  cal¬ 
culated  by  Eqn .  (22).  This  calculation  ignores  inertial 

effects  and  includes  both  static  and  dynamic  deformation 


Flexural  Response  of  a  90°  Unidirectional  APC-2  Specimen 


A  similar  response  is  shown  in  Figure  5  for  the  neat 
resin  specimen.  This  response  correlates  well  with  the 
nearly  t empe ra t ur e - i nde pe nde n t  linear  tensile  response  ob¬ 
served  by  Cebo  for  PEEK  films  below  strains  of  1%  (4:495). 

The  regions  of  the  P  -  wQ  curves  used  for  dynamic 
displacements  are  also  outlined  in  both  figures.  In  neither 
case  did  the  total  of  static  and  dynamic  displacements 
exceed  the  linear  region.  Extensive  pretest  flexure  tests 
determined  the  best  S/N  ratio  for  the  90°  beam.  This  opti¬ 
mum  value  was  obtained  with  a  static  to  dynamic  displacement 
ratio  of  about  five  to  one.  The  static  displacement  applied 
to  90°  samples  was  0.25mm,  which  resulted  in  a  low-range 
displacement  transducer  output  of  5  volts.  This  is  the 
center  of  the  dynamic  range  depicted  in  Figure  4.  Similar 
tests  were  conducted  for  the  0°  beam  and  neat  resin  samples, 
resulting  in  static  to  dynamic  displacement  ratios  of  five 
to  three  for  each  . 

Initially,  the  dynamic  mechanical  behavior  of  the  two 
unidirectional  composites  will  be  presented  here  separately 
for  the  sake  of  clarity.  The  first  to  be  discussed  will  be 
the  matrix-dominated  geometry,  i.e.,  the  90°  specimen.  How¬ 
ever,  a  word  on  notation  is  first  required. 

Note  that  in  Figure  6,  and  all  other  plots  generated  by 
the  Dynastat,  the  in-phase,  or  storage  modulus,  is  labeled 
as  M',  not  E'.  The  Dynastat  reserves  the  symbols  E'  and 
E' '  for  unsealed  stiffness.  The  units  of  M'  and  M''  are 


identical  to  those  of  E'  and  E'*  from  Eqn.  (20).  Of  course, 
the  dimensionless  loss  factor,  tan  j  ,  or  E"/E',  is 

the  same  regardless  of  notation. 

Trans  verse  (  9  0  °  J_  Behavior 

The  90  sample  underwent  three-point  flexure  from  2 3 °C 
to  280°C.  The  temperature  response  at  1  rad/sec  is  depicted 
in  Figure  6.  For  the  same  specimen  and  test,  the  responses 
at  the  frequencies  of  10  rad/sec  and  100  rad/sec  are  pre¬ 
sented  in  Figures  7  and  8.  The  following  descr i pt i ons  apply 
to  all  three  frequencies. 

The  in-phase  component  is  seen  to  remain  fairly  con¬ 
stant  with  temperature  until  120° C.  A  /}  transition  may 
occur  between  40°C  and  7  0  °C .  More  likely,  this  is  noise. 

In  the  glassy  region,  the  in-phase  component  is  very  large, 
and  the  loss  modulus  is  extremely  small.  Because  of  this, 
tan  (J  approaches  the  machine  lower  measurement  limit.  Data 
in  this  region  is  thus  more  susceptible  to  noise  and  the 
scatter  in  the  small  loss  modulus  values  represents  noise 
more  than  a  material  effect. 

The  loss  modulus  values  reach  a  broad  and  shallow  peak 
in  the  transition  region.  The  out-of-phaso  response  beyond 
the  transition  region  declines,  and  then  increases  as  the 
temperature  nears  the  crystalline  melt  temperature  of  30F°c 
observed  by  Tung  and  Dynes  (5:508)  .  In  this  rubbery  region, 
the  storage  modulus  is  seen  to  steadily  decline. 
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Temperature  Dependency  of  the  Complex  Flexural  Moduli  of  APC-2  90°  Specimen 


Shift  factors  for  the  90  specimen  are  displayed  in 
Figure  9,  relative  to  the  reference  temperature  of  210. 4°C. 
The  shift  factors  will  vary  with  reference  temperature,  of 
course,  but  Figure  9  indicates  that  they  also  depend  on 
which  variable  ( M ’  or  M" )  is  shifted.  Those  calculated  by 
the  Datalyzer  are  labeled  "Super"  because  that  is  the  soft¬ 
ware  command  'equired  by  the  Datalyzer  to  begin  shift  factor 
calculation-*. 

Figure  9  shows  shift  factors  calculated  by  other  means, 
too.  Those  labeled  "WLF"  were  calculated  from  Eqn .  (21). 

Those  shift  factors  labeled  "Hand"  were  determined  manually 
by  plotting  individual  spectra  and  passing  them  across  a 
scaled  light  table.  It  is  estimated  that  the  precision  of 
the  manual  shifts  due  to  the  grid  size  is  no  better  than  3%. 
Due  to  inexperience  and  the  difficulty  with  visually  mini¬ 
mizing  the  area  between  two  overlapping  curves,  the  typical 
accuracy  is  close  to  20%. 

Because  of  the  inaccuracy  of  the  manual  method,  and  due 
to  the  possible  poor  fit  of  the  polynomial  spline  at  the 
edges  of  individual  spectra,  trends,  and  not  values,  should 
be  compared  in  Figure  9.  The  WLF  equation  was  applied  with 
"universal"  constants.  No  attempt  was  made  to  calculate  new 
constants,  unique  to  APC-2,  in  the  WLF  equation  from  slope 
and  intercept  values. 

The  general  trend  of  the  experimental  shift  factors, 
corresponding  to  the  storage  modulus,  whether  curve-fit  by 


Fig.  9.  Effect  of  Reduction  Technique  on  90 


hand  or  by  the  computational  method  of  the  Datalyser,  chows 
agreement  with  a  t he r mor he o log i ca 1 1 y  simple  solid  that  the 
WLF  equation  models.  The  leveling  of  the  log ^  a(T)  at 
lower  temperatures  indicates  that  the  WLF  equation,  as  fit 
to  this  data,  should  be  modified  to  include  an  energy  term, 
similar  to  an  Arrhenius  expression  for  the  activation  energy 
(11:290). 

The  factors  based  on  shifts  of  the  loss  modulus  are 
small,  and  "Super"  has  difficulty  detecting  any  trend.  This 
is  thought  to  be  due  to  the  noisy  and  relatively  featureless 
loss  modulus  spectra.  The  cur ve - f i tt i ng  algorithm  of  the 
Fafalyser  is  probably  too  sensitive  to  large  changes  in 
amplitude  or  to  large  differences  in  the  slopes  at  the  ends 
of  the  splines.  Hand  shifting  was  extremely  difficult. 

O 

Long  i  tud  i  na  1  (  0  Be  ha  vi  or 

Returning  to  the  discussion  of  the  fiber  dominated,  or 
0°,  specimen.  Figure  10  displays  the  temperature  dependency 
of  the  composite  over  the  range  from  80° C  to  230°C.  Two  at¬ 
tempts  to  obtain  data  at  higher  temperatures  were  thwarted, 
once  due  to  mechanical  problems,  and  once  to  a  loss  of  RAM 
data.  Note  the  difference  in  the  vertical  scaling  between 
the  loss  and  storage  moduli,  chosen  for  convenience. 

As  expected,  the  storage  modulus  is  extremely  high,  two 
orders  of  magnitude  above  the  out-of  phase  component,  and 
shows  little  temperature  or  frequency  dependence.  A  transi¬ 
tion  region  for  the  storage  modulus  is  barely  apparent  on 
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the  scale  of  Figure  10.  Over  the  entire  temperature  range, 
the  storage  modulus  decreases  by  less  than  4%. 

The  most  striking  feature  of  the  fiber  dominated  behav¬ 
ior,  however,  is  the  depressed  region  around  T  for  the  loss 

§ 

modulus.  Only  the  smallest  peaks  above  this  "valley"  are 
discernible  at  either  frequency  shown  in  Figure  10.  The 
trend  above  this  transition  range,  though,  matches  the  be¬ 
havior  of  the  90  °  sample  in  that  the  loss  factor  continues 
to  rise  with  temperature.  Significantly,  the  loss  factor  is 

O 

typically  an  order  of  magnitude  less  than  that  of  the  °0 
sample . 

O 

The  shift  factors  generated  from  the  spectra  of  the  0 
specimen  are  displayed  in  Figure  11  as  a  function  of  temper- 

O 

iture,  again  with  respect  to  a  reference  of  210  C.  The 
polynomial  fit  of  "Super"  exaggerates  the  shift  factor  for 
the  storage  modulus  at  several  points,  but  matches  the  WLF 
curve  more  closely  than  the  hand  fit.  Hand -sh  i  f  t  i  r.g  was 


extremely  difficult  due  to  the  small  changes  in  slopes  and 
magnitudes  of  individual  spectra,  but  yielded  a  master  curve 
with  better  continuity  over  a  larger  shifted  spectrum.  The 
master  curve  will  be  discussed  below.  As  with  the  90°  beam, 
shift  factors  based  on  the  loss  modulus  are  small,  and  their 
trend  is  flat,  regardless  of  the  technique  employed. 


e  r 


Master  Curves 


Several  master  curves  were  drawn  for  the  two  specimen 
fiber  orientations.  In  each  case,  the  shift  factors  used 
are  identified  by  computational  method:  by  hand,  or  by 
eighth-order  polynomial  spline  ("Super").  The  moduli  are 
weighted  by  the  ratio  of  the  reference  temperature  to  actual 
temperature  at  which  the  data  was  obtained. 

Figure  12  is  the  master  curve  for  the  90°  specimen 
generated  by  the  Datalyzer  "Super"  routine  for  the  storage 
modulus.  Immediately  noticeable  are  several  gups  in  the 
spectrum  at  the  reduced  frequencies  above  100  GHz.  These 
gaps  in  turn  produce  a  huge  overall  bandwidth  of  the  reduced 
frequency  spectrum.  This  was  caused  by  carrying  out  the 
frequency -temperature  superposition  for  all  temperatures., 
despite  the  assertion  at  the  beginning  of  this  chapter.  The 
results  clearly  indicate  the  fallacy  of  extending  superposi¬ 
tion  into  the  glassy  region  with  noisy  data. 

Even  in  the  applicable  temperature  regime,  however,  the 
large  shift  factor  near  the  reference  temperature  leaves  a 
gap  in  the  master  curve  of  the  storage  modulus  just  below  1 
GHz.  While  the  master  curves  at  lower  temperatures  should 
not  be  displayed,  this  particular  gap  indicates  that  "Super" 
must  still  be  supplemented  with  manual  entries  by  the  inves 
t i ga t  or  . 

Tbs'  remaining  portions  of  the  master  curves  in  Figure 
12  resemble  the  complete  curves  presented  in  Figure  13,  the 


master  curves  generated  by  hand  for  the  storage  modulus. 

Figure  13,  shifting  was  applied  only  near  and  above  T  . 

8 

storage  modulus  can  be  seen  to  decline  at  the  lower  frequ 
cies  in  a  manner  similar  to  other  highly  crystalline  poly 
me  r  s  (11:42). 

The  loss  modulus  is  very  low  in  the  region  where  the 
storage  modulus  changes  slowly  with  frequency.  The  wide 
scatter  at  higher  frequencies  may  be  due  to  the  presence 
the  fibers,  but  it  also  confirms  the  assertion  that  one 
can't  apply  frequency-temperature  superposition  to  crysta 
line  polymers  in  the  glassy  state  (29:199).  The  broad 
transition  region  correlates  with  the  observations  at  low 
strains  of  other  polymers  with  an  intermediate  degree  of 
crystallinity  (29:219).  The  loss  modulus  increases  with 
decreasing  frequencies  following  a  decline  below  the 
plateau-like  transition  region.  This  trend  indicate::  the 
increasing  viscoelastic  nature  of  the  composite  above  T  . 

O 

The  master  curves  for  the  90  specimen  developed  by 
generating  shifts  relative  to  the  loss  moduli  spectra  are 
presented  in  Figures  14  and  15.  The  folding  effect,  or 
reversal  of  the  master  curve  at  low  frequencies,  from 
"Super"  shifts  is  evident  in  Figure  14.  This  curve  has 
little  value. 

In  Figure  15,  the  loss  modulus  was  used  as  a  refer  on 
f'r  hand  shifting  spectra  only  at  temperatures  above  i  ",p° 
Again,  the  broad,  shallow  peak  in  the  loss  modulus  appear 
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:  n  the  transition  region.  The  loss  modulus  decreases  wit 
lower  frequencies  from  this  peak,  and  then  increases  wifi- 
higher  temperatures/lower  frequencies. 

The  behavior  of  both  moduli,  and  the  loss  factor, 
resembles  the  master  curves  constructed  in  Figure  13.  F  i 
foronces  are  seen  mainly  in  the  extension  of  the  scale  cf 
the  reduced  frequency.  The  master  curves  generated  from 
shift  factors  based  on  the  storage  modulus  cover  a  wider 
range.  f  frequencies  .  The  master  curves  from  the  loss 
modulus  shift  factors  indicate  the  range  of  values,  or 
.-Titter,  at  a  particular  frequency. 

The  same  approach  as  was  just  presented  in  detail  fc 

0  o 

the  00  specimen  produced  master  curves  for  the  0  specie 
Figure-  16  and  17  display  the  results  of  shifting  the  a 1 0 
age  modulus,  with  respect  to  a  reference  temper  a fu re  of 

O 

210  C,  by  "Super"  and  by  hand.  As  with  the  master  curves 
the  ?°°  specimen,  "Super"  produces  gaps  of  several  decade 
and  folds  some  spectra  over  others.  There  is  1  i  1 1 1  a  cor 
relation  of  the  resulting  master  curve  to  any  trend.  Th.  i 
demonstrates  the  inability  of  the  Datalyser  to  handle  no: 
or  featureless  data. 

Shifting  the  data  by  hand,  however,  producer  the  lr ' 
ohall-W  peak  in  the  cut  -of  -phase  response  s--'en  previous] 
‘-he  no  0  .-simple.  The  range  of  values  that  the  loss  rr.edulu 
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limit  of  the  capability  Dynastat  to  resolve  mechanical  re¬ 
sponse  from  noise.  This  behavior  is  expected  for  the  fiber- 

o 

dominated  0  specimen.  The  storage  modulus  declines  only 
gradually  through  the  transition  region,  and  is  nearly  con¬ 
stant  over  a  wide  spectrum  of  low  frequencies. 

When  the  data  is  shifted  with  respect  to  the  loss 
modulus,  "Super"  again  folded  the  master  curve  over  itself 
near  the  reference  spectra  between  1  Hz  and  10Hz.  Figure  18 
shows  this  behavior,  and  is  of  little  value. 

The  manual,  or  hand,  shift  factors  produced  the  master 
curve  presented  in  Figure  19.  This  curve  shows  the  compres¬ 
sion  of  the  reduced  frequency  scale  due  to  the  difficulty  of 
shifting  noisy  data  by  hand.  It  too  is  of  little  value. 

Shear  Deformation 

Equation  15a  predicts  the  flexural  modulus  with  the 
contribution  of  shear  deformation  included.  Table  3  incor¬ 
porates  the  values  for  the  shear  modulus  of  APC-2  determined 
at  80°C  and  at  190°C  by  Tung  and  Dynes  from  Ref.  5.  The 
experimentally  determined  values  for  the  0°  specimen  of  the 
effective  bending  modulus,  or  complex  modulus,  at  the  same 
frequency  as  used  in  Ref.  5,  are  also  listed.  These  values 
are  then  increased  by  the  method  described  in  the  second 
chapter  of  this  report.  The  results,  tabulated  in  the  far 
right  column,  approximate  the  parts  of  the  effective  flex¬ 
ural  modulus  due  to  bending  alone  for  the  given  load.  Shear 
deformation  is  seen  to  contribute  to  the  total  modulus  by 
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less  than  8%  at  or  below  190  C  for  these  typical  examples. 
Note  Eqn.  (15a)  predicts  an  increase  in  E  with  temperature. 
This  isn't  the  general  case.  Rather,  G  is  so  small  around 
190JC  that  Eqn.  (15a)  understates  the  contribution  of  shear 
deformation  to  the  flexural  modulus. 


Table  3 

Comparison  of  the  Experimental  Complex  Flexural  Modulus 
With  Shear  Deformation  Theory 


Temperature 

(°C) 


Maximum 
Center 
Span  Load 
(N) 


Experimental 

Flexural 

Modulus 

(GPa) 


Shear 

Modulus 

(GPa) 


Predicted 

Flexural 

Modulus 

(GPa) 


80 


7.57  136.8 


4.9  139.6 


190 


7.38  133.3 


2.28  143.5 


Comparison  to  Rule  of  Mixtures  Approach 

The  abbreviated  frequency-temperature  sweep  conducted 
on  the  PEEK  neat  resin  flexure  sample  obtained  slightly 
ambiguous  results.  High  temperature  data  is  extremely  sus¬ 
pect.  The  output  of  the  load  cell  appeared  as  a  clipped 
trace  on  the  oscilloscope  during  some,  but  not  all,  of  the 
cycles  at  each  frequency  at  210. 8°C.  The  indentor  rose  off 
the  sample  at  these  points.  This  behavior  was  not  repeated 
in  later  tests  on  composite  samples.  The  difficulty  in 
producing  neat  resin  samples  prevented  a  repeat  in  this 


test . 


It  is  not  possible  to  say  whether  the  cause  of  the 
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the  indcntor . 


Confidence  in  the  lower  tempera  tor' 


:  h  o  u '  j  h 


Figure  20  shows  the  thermal  dependence  of  t!> 


tempi ox  flexural  moduli  of  the  PEEK  neat  resin  beam.  M 


ml  below  T  ,  ther< 
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spread  b  e  t ■ 
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".emper  ature  .  The  overall  behavior  of  t!v.  material 
that  of  the  matrix-dominated  geometry  (00°  sample).  Tl: 
; ts  s  y  region  is  evident,  and  a  1  ar  go  pea!:  ocrur s  :  r.  th 
out f -phase  modulus  near  1^0°C.  The  decrease  in  the  1 
modulus  above  T  is  less  than  the  increase  from  the  gin 
:::,  as  also  exhibited  by  the  composite. 

The  complex  moduli  cor  re  spending  to  some  of  the  dr 
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and  for  the  length-to-vidth  ratio  of  L/b  =  3.33,  the  maximum 
corrections  become  factors  of  1.035  for  an  isotropic  beam 
and  1.91  for  unidirectional  APC-2. 

Thus,  by  this  method,  the  predicted  value  of  the  effec¬ 
tive  flexural  modulus  of  APC-2  would  differ  from  the  tensile 
modulus  by  almost  a  factor  of  two.  This  occurs  where  i>  is 
large,  i.e.,  at  higher  temperatures.  The  discrepancy  also 
shows  the  importance  of  properly  sizing  the  flexure  speci¬ 
men,  as  noted  previously  (Reference  6).  Mair  concluded  that 
the  discrepancy  decreased  with  increasing  L/b  ratios.  The 
Dynastat  test  chamber  was  restricted  in  span,  as  mentioned 
in  the  previous  chapter.  Only  by  reducing  the  width  would 
the  effects  of  shear  coupling  become  less  noticeable,  but 
this  would  have  introduced  the  possibility  of  increased  free 
edge  effects. 

Thus,  shear  coupling  was  not  completely  removed  from 
the  beam  response.  Experimental  values  otherwise  show  good 
agreement.  One  should  observe,  though,  that  Eqn.  7  should 
not  be  applied  to  angle-ply  laminates,  and,  indeed,  even 
corrections  like  Eqn.  (7)  may  overestimate  the  flexural 
modulus  by  a  factor  of  six  (8:40). 
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S umma r  v  of  Resul  ts 

An  analysis  of  three-point  flexure  for  the  viscoelastic 
character i zat i on  of  a  laminated  composite  of  polyether ether  - 
ketone  (PEEK)  has  been  conducted  based  on  an  extension  of 


laminated  plate  theory  and  the  work  of  Whitney,  Browning, 
and  Mair.  The  results  of  a  dynamic  mechanical  analysis  of 
the  thermoplastic  matrix/graphite  fiber  composite  APC-2  have 
been  compared  to  the  analytic  model.  Specimen  geometry  was 
selected  to  avoid  any  considerable  warping  due  to  shear 
deformation  and  biaxial  response  due  to  shear  coupling. 

Flexural  data  was  obtained  with  a  Dynastat  viscoelastic 
dynamic  and  static  analyzer  over  a  two-decade  frequency 
range  and  from  roughly  ambient  temperatures  to  more  than 
230  C.  Low  strain  levels  produced  a  macroscopic  linear 
stress  behavior.  Experimental  flexural  moduli  at  tempera¬ 


tures  below  190°c,  reduced  from  a  model  that  neglects  shear, 
match  within  5%  the  predictions  that  include  shear  deforma¬ 
tion.  The  T„  was  measured  to  be  between  150°C  and  160°C  for 
frequencies  from  0.16  Hz  (1  rad/sec)  to  15.90  Hz  (100 
rad/sec ) . 

Master  curves  have  been  created  for  the  frequency- 
temperature  behavior  of  two  unidirectional  geometries:  one 
with  the  reinforcing  fibers  parallel  with  the  longitudinal 
axis,  the  other  with  the  fibers  oriented  transversely.  The 
shift  factors  which  generated  the  curves  have  been  compared 


by  technique  and  judged  for  their  effectiveness.  The  shift 
factors  computed  by  the  Dynastat  for  the  storage  modulus 
compare  reasonably  well  with  the  WLF  equation  with  standard 
constants.  The  Dynastat  produced  master  curves  similar  to 
those  created  from  shift  factors  determined  by  hand  shifting 
the  experimental  data. 

The  effect  of  matrix  viscoelasticity  on  the  overall 
composite  behavior  was  demonstrated  in  the  master  curves. 

The  composite  demonstrated  nearly  constant,  f iber -dominated 
flexural  rigidity  even  in  the  transverse  geometry,  but  be¬ 
cause  of  frequency-  and  temperature-dependent  changes  in  the 
loss  factor  of  one  magnitude  in  order,  the  damping  charac¬ 
teristics  of  the  composite  vary  significantly.  Broad,  shal¬ 
low  peaks  in  the  loss  modulus  and  the  loss  factor  in  a 
transition  zone  correlate  well  with  neat  resin  flexural 
behavior.  After  declining  with  lower  frequencies  from  the 
transition  zone,  the  out-of-phase  components  increased  at 
extremely  low  frequencies  in  a  manner  indicating  the  vis¬ 
coelastic  effect  increased  with  temperature  above  T  .  These 
trends  were  observed  in  both  geometries  at  roughly  the  same 
frequencies  and  tempera tur es . 

Neat  resin  response  in  flexure  at  a  smaller  number  of 
temperatures  and  frequencies  was  obtained  and  incorporated 
in  a  rule  of  mixtures  approach  to  predict  composite  behav¬ 
ior.  Where  high  quality  neat  resin  data  was  available,  the 
predictions  compared  well  with  the  experimental  results. 


Conclusions 


The  objectives  of  this  research  are  repeated  here. 

Objective  1.  Evaluate  the  flexural  test  and  account 
for  shear  and  Poisson  effects.  The  viscoelastic  behavior  of 
APC-2  is  predicted  adequately  by  adapting  elastic  static 
plate  theory  to  the  dynamic  flexure  of  unidirectional  lami¬ 
nates.  The  degree  of  APC-2  anisotropy  is  not  significant  a t 
the  load  levels  within  this  research  .  Shear  deformation, 
Poisson  effects,  and  shear  coupling  do  not  contribute  sig¬ 
nificantly  to  the  effective  flexural  modulus  determined  by 
the  three  point  bend  test  if_  one  properly  sizes  the  specimen 
geometry  and  uses  small  de  format  ions  as  i n  this  research  . 

Object  i  ve  2.-  Characterize  the  frequency  and  tempera¬ 
ture  dependency  of  APC-2  in  dynamic  flexure.  The  complex 
flexural  modulus  of  APC-2  reflects  a  strong  fiber  dominance 
which  suppresses  the  out-of-phase  component.  Master  curves 
of  un id irectional  laminates  are  extremely  difficult  to  ob¬ 
tain.  Best  results  are  obtained  for  the  matrix  dominated 
flexure  geometry  (90°  sample).  When  based  on  the  behavior 
of  the  storage  modulus,  the  master  curves  do  exhibit  the 
expected  viscoelastic  nature  of  the  thermoplast ic  matrix. 

Master  curves  should  not  be  generated  from  the  behavior 
of  the  loss  modulus  of  composite  materials.  The  loss  modu¬ 
lus  data  must  be  of  high  quality  and  exhibit  marked  fre¬ 
quency  and  temperature  characteristics. 


The  presence  of  fibers  in  the  matrix  dominated  orienta¬ 
tion  slightly  increases  the  glass  transition  temperature  of 
APC-2  from  the  Tg  of  PEEK.  Otherwise,  the  effective 
flexural  modulus  of  this  thermoplastic  composite  can  be 
predicted,  from  a  rule  of  mixtures  approach,  if  quasi-static 
behavior  is  assumed,  the  effects  of  crystallinity  are  con¬ 
sidered,  and  the  unidirectional  fibers  are  oriented  trans¬ 
verse  to  the  neutral  axis  of  the  beam. 

The  presence  of  fibers  alters  the  temperature  depend¬ 
ency  of  the  out-of-phase  component  of  the  complex  flexure 
modulus  for  the  composite  with  fibers  aligned  parallel  with 
the  neutral  axis.  The  viscoelastic  effect  is  unexpectedly 
high  in  the  glassy  region.  This  behavior  is  currently 
unexplained,  and  merits  further  research. 

Object i ve  2*  Compare  the  f r equency- temperature  shift 
factors  obtained  in  three  ways:  analyt ically,  using  the  WLF 
equation;  experimentally,  using  the  Dynastat  cur ve - f i tt i ng 
algorithm;  and  manually,  using  hand-shifted  values  from 
experimental  results.  As  mentioned  previously,  the  highest 
quality  master  curves  were  produced  by  the  Dynastat  by 

O 

shifting  the  storage  modulus  spectra  for  the  90  sample. 
These  shift  factors  most  closely  resemble  the  behavior  of  a 
thermorheologically  simple  solid  that  the  WLF  equation 
models.  The  difference  between  the  Dynastat  shift  factors 
and  those  obtained  manually  can  be  expected  to  decrease  with 
operator  experience  with  the  manual  technique.  However,  the 


Dynastat  can  provide  reliable  and  good  resolution  data.  The 
automatic  curve-fitting  capability  is  useful  as  a  first  step 
in  calculating  f r equency- temperatur e  superpositions.  Its 
utility  in  reducing  the  time  required  to  generate  shift 
factors,  and  plot  the  subsequent  master  curves,  is  evident. 
The  results,  however,  must  be  obtained  from  test  conditions 
known  to  produce  linear  behavior  with  a  high  signal  to  noise 
ratio. 

Recommendations 

The  analysis  employed  in  this  study  appears  valid  for 
unidirectional,  anisotropic,  laminated  composite  materials 
with  high  flexural  rigidities.  However,  future  studies 
should  employ  some  experimental  modifications: 

1.  Flexure  tests  should  cover  a  wider  range  of  temper¬ 
atures  than  used  in  this  research.  The  cryogenic  and  near¬ 
melt  temperature  regions  were  not  included  here,  and  much 
useful  design  information  can  be  obtained  from  their  inclu¬ 
sion  into  a  test  plan. 

2.  Apart  from  an  increase  in  the  signal  to  noise  ratio 
that  the  experimental  method  produces,  frequency- temperature 
shifting  requires  data  with  better  features.  Combinations 
of  the  following  steps  may  improve  the  capability  of  the 
Dynastat  and  the  three  point  flexure  test  to  obtain  visco¬ 
elastic  characteristics  of  composite  laminates.  The  temper¬ 
ature  increment  should  be  decreased  near  Tg  ,  or  the  number 
of  discrete  frequencies  used  at  each  temperature  should  be 


increased.  This  should  yield  polynomial  curve  fits  of  the 
modulus  spectra  that  will  be  less  susceptible  to  noise  and 
slope  changes  at  the  ends  of  the  splines. 
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3.  The  Materials  Laboratory  should  investigate  the 
feasibility  of  using  four  point  flexure  in  future  viscoelas¬ 
tic  character izations .  This  loading  subjects  the  center 
span  of  the  flexure  specimen  to  a  constant  bending  moment 
and  more  accurately  resembles  a  state  of  pure  bending.  The 
three  dimensional  state  of  stress  at  the  supports  and  inden- 
tors  will  of  course  still  exist,  but  it  may  not  introduce  as 
many  complications  to  the  analysis. 

4.  Tensile  tests  of  APC-2  coupons  may  better  define 
the  frequency-temperature  dependence  of  the  Poisson  ratio. 
Shear  coupling  effects  of  other  laminates  to  be  tested  on 
the  Dynastat  in  the  future  may  be  better  controlled  by  this 
knowledge . 


5.  The  Materials  Laboratory  should  investigate  the 
feasibility  of  using  test  chamber  supports  with  longer  spans 
for  the  Dynastat.  This  would  additionally  minimize  shear 
coupling  and  free  edge  effects. 

6.  Finally,  the  test  technique  to  be  used  with  the 
Dynastat  must  improve  the  signal  to  noise  ratio  observed  in 
this  research.  Observations  over  15  cycles  should  be  used 
with  stiff  composites,  while  a  smaller  number  of  cycles  may 
still  be  adequate  with  homogeneous  isotropic  viscoelastic 
materials.  If  necessary,  higher  deformations  should  be  used 
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compared  to  shift  factors  predicted  from  the  Williams- 
i.ande  1-Ferrv  equation.  Determining  composite  behavior  from 
knowledge  of  matrix  and  fiber  behavior  is  still  applicable. 
The  capability  of  the  Dvnastat  to  perform  viscoelastic 
characterizations  of  stiff,  unidirectional  composites,  base 
on  comparisons  of  the  various  shift  factors,  is  reduced  due 
to  the  presence  of  the  reinforcing  fibers.  Recommendations 
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ilvnamic  s  i  g n a  1  -  t  o - n  o  i  se  ratio,  are  described. 
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